The first fermi in a high energy nuclear collision. 
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OAt very high energies, weak coupling, non— perturbative methods can be used to 
study classical gluon production in nuclear collisions. One observes in numerical 
f*\ ' simulations that after an initial "formation" time, the produced partons are on 

shell, and their subsequent evolution can be studied using transport theory. At 
the initial formation time, a simple non— perturbative relation exists between the 
energy and number densities of the produced partons, and a scale determined by 
K ^ , the saturated parton density in the nucleus. 

*yv . An outstanding problem in high energy scattering is the problem of ini- 

tial conditions for particle production El. In perturbative QCD, for processes 

t— I . which involve a hard scale Q 2 ^> A.q CD , the hard and soft contributions can 

be factorized. The soft contributions are lumped into non-perturbative, pro- 
cess independent parton distribution functions, while the hard contributions 

,S3 ' are computed for each physical process of interest. For a fixed hard scale of 

interest Q 2 , there is a center of mass energy yfs beyond which this approach 

Q, 1 in particular, and the operator product expansion (OPE) in general breaks 

down a. However, since the parton densities, in this regime are large, weak 

i-^ ' coupling classical methods may be applicable 0. Wilson renormalization group 

^- . methods have been developed for this high parton density regime □. 

At small cc, classical parton distributions in a nucleus can be computed 
in a model with a dimensionful scale /i 2 proportional to the gluon density per 
unit transverse area. In this model, parton distributions saturate at a scale 
Q s oc g 2 /i. For Au-Au collisions at RHIC, one can estimate Q s ~ 1 GeV, 
and at the LHC, the saturation scale will be Q s ~ 2-3 GeV. Most of the 
gluons produced therefore have transverse momenta k t ~ Q s , and since this 
scale for RHIC and LHC is at least marginally a weak coupling scale, these 
classical methods may be applied to study the production and initial evolution 
of partons at RHIC and LHC. 

These classical methods were first applied to nuclear collisions by Kovner, 
McLerran and WeigertH. For an interesting alternative approach, see Ref. U. 
Assuming boost invariance, and matching the equations of motion in the 
forward and backward light cone, they obtained the following initial condi- 
tions for the gauge fields in the A T = gauge: A l ± \ T= a — A\ + A\, and 
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A ± \ T=0 = ± l -^x ± [A{,A^}. Here A\ 2 (p ± ) (i = 1,2) are the pure gauge trans- 
verse gauge fields corresponding to small x modes of incoming nuclei (with 
light cone sources p^ S(x T )) in the 9(±x~)8(^fx + ) regions respectively of the 
light cone. 

The sum of two pure gauges in QCD is not a pure gauge-the initial condi- 
tions therefore give rise to classical gluon radiation in the forward light cone. 
For pt » asp, the Yang-Mills equations may be solved perturbatively to 
quadratic order in asp/pt- After averaging over the Gaussian random sources 
of color charge p ± on the light cone, the perturbative energy and-, number 
distributions of physical gluons were computed by several authors uu. In the 
small x limit, it was shown that the classical Yang-Mills result agreed with 
the quantum Bremsstrahlung result of Gunion and BertschB. 

In Ref.Q, we suggested a lattice discretization of the classical EFT, suitable 
for a non-perturbative numerical solution. Assuming boost invariance, we 
showed that in A T = gauge, the real time evolution of the small x gauge fields 
A±(xt,r), A n (xt,T) is described by the Kogut-Susskind Hamiltonian in 2+1 
dimensions coupled to an adjoint scalar field. The lattice equations of motion 
for the fields are then determined straightforwardly by computing the Poisson 
brackets. The initial conditions for the evolution are provided by the lattice 
analogue of the continuum relations discussed earlier in the text. We impose 
periodic boundary conditions on an N x N transverse lattice, where TV denotes 
the number of sites. The physical linear sizajpf the system is L — a N, where a 
is the lattice spacing. It was shown in Ref.El that numerical computations on 
a transverse lattice agreed with lattice perturbation theory at large transverse 
momentum. For details of the numerical procedure, and other details, we refer 
the reader to Ref.Ej. 

In our numerical simulations, all the relevant physical information is com- 
pressed in g 2 p, and L, and in their dimensionless product g 2 pLtA. The strong 
coupling constant g depends on the hard scale of interest; p ex A 1 / 6 depends 
on the nuclear size, the center of mass energy, and the hard scale of interest; 
L 2 is the transverse area of the nucleus. Assuming g = 2 (or as = 1/tt), 
H = 0.5 GeV (1.0 GcV) for RHIC (LHC), and L = 11.6 fm for Aw-nuclei, 
we find g 2 pL sa 120 for RHIC and sa 240 for LHC. (The latter number would 
be smaller for a smaller value of g at the typical LHC momentum scale.) As 
will be discussed later, these values of g 2 pL correspond to a region in which 
one expects large non-perturbative contributions from a sum to all orders in 
~ 6asp/pt, even if as <C 1. We should mention here that deviations from 
lattice perturbatkm theory, as a function of increasing g 2 pL, were observed in 
our earlier workEj. 

In Ref. t3, we computed the energy density e as a function of the proper 



time r. This computation on the lattice is straightforward. To obtain this re- 
sult, we computed the Hamiltonian density on the lattice for each p , and then 
took the Gaussian average (with the weight p 2 ) over between 40 p trajectories 
for the larger lattices and 160 p trajectories for the smallest ones. 

The dependence of et as a function of r was investigated in our numer- 
ical simulations. For larger values of g 2 pL, et increases rapidly, develops a 
transient peak at r ~ l/g 2 p, and decays exponentially there onwards, satis- 
fying the relation a + [3e~' yT , to the asymptotic value a (equal to the lattice 
dE/L 2 /drjl). This behavior is satisfied for all g 2 pL > 8.84, independently of 
N. One can interpret the decay time td — Vt/s 2 ^ as the appropriate scale 
controlling the formation of gluons with a physically well defined energy In 
other words, tjj is the "formation time" in the sense used by Bjorkenliil 
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Figure 1: £t/(<ji 2 /j) 3 extrapolated to the continuum limit: / as a function of g 2 fJ,L. The 
error bars are smaller than the plotting symbols. 



The physical energy per unit area per unit rapidity of produced gluons can 
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be defined in terms of a function f(g 2 fiL) as 

1 dE _ 1 
L 2 drj g 2 



f(g 2 nL) {g 2 tf • (1) 



The function / here is obtained by extrapolating our results for finite lattice 
spacings to the continuum limit. In the region of physical interest for heavy 
ion collisions, / varies very slowly. It changes by ~ 25% for nearly an order of 
magnitude change in g 2 iiL. The saturation scale Q s ~ 6 as/x-one can therefore 
re-write our result for the energy density in terms of Q s . 

Doing so, wejjcpnfirmed that our results are consistent with an estimate 
by A. H. Mueller Ej for the number of produced gluons per unit area per unit 
rapidity. He obtains dN/L 2 /drj = c (N 2 ~l) Ql/4n 2 as N c , and argues that the 
number c is a non-perturbative constant of order unity. If most of the gluons 
have p t ~ Q s , then dE/L 2 jdr\ = c' (N 2 — 1) Q^/4n 2 as N c which is of the same 
form as our Eq. [l[ In the g 2 [iL region of interest, our function / rs 0.23-0.26. 
We obtain c' = 4.3-4.9. Since one expects a distribution in momenta about Q s , 
it is very likely that c' is at least a factor of 2 greater than c-thereby yielding 
a number of order unity for c as estimated by Mueller. This coefficient can 
be determined more precisely when we compute the non-perturbative number 
and energy ^distributions £3. 

In Ref. £3, we estimated the initial energy per unit rapidity of produced 
gluons at RHIC and LHC energies. We did so by extrapolating from our SU(2) 
results to SU(3) assuming the N c dependence to be (N 2 — l)/N c as in Mueller's 
formula. At late times, the energy density is e = {g 2 ^) A f(g 2 nL)-f(g 2 [iL)/g 2 , 
where the formation time is td = l/ r y(g 2 ^iL)/g 2 fj, as discussed earlier. Wc 
find e RHIC « 66.49 GeV/fm 3 and e LHC « 1315.56 GeV/fm 3 . Multiplying 
these numbers by the initial volumes at the formation time td, we obtained 
the classical Yang-Mills estimate for the initial energies per unit rapidity Et 
to be E$ HIC « 2703 GeV and E^ HC « 24572 GeV respectively. 

Compare these numbers to results presented recently by KajantieE3 for the 
mini-jet energy (computed for p t > p sa tr-wheTe p sa t is a saturation scale akin 
to Q s ) in the pQCD mini-jet approachM He obtains E^ HIC = 2500 GeV 
and E^ HC = 12000. The remarkable closeness between our results for RHIC is 
very likely a coincidence. Kajantie's result includes a K factor of 1.5-estimates 
range from 1.5-2.5 E£l For the. latest estimates from our Finnish colleagues, see 
the preprint of Eskola et al.E3. If we pick a recent value of K rs 2E3, we obtain 
as our final estimate, E$ HIC w 5406 GeV and E^ HC w 49144 GeV. 

We can also boldly estimate the number of produced gluons at central 
rapidities. As mentioned in the preceding text, the value of the constant c 
in the expression for the number distribution is currently being computed 



numerically. We obtain for Au-Au collisions in one unit of rapidity the result 
that NrhiC — 714 • c and Nlhc = 2855 • c. Given that the corresponding 
constant for the energy density was larger, we would anticipate that it is more 
likely that c = 2-2.5. Taking the higher value, we obtain Nrhic = 1785 and 
Mlhc — 7138. Again, these values are very close to those of Eskola et al. 
E3. Note that they consider Pb-Pb collisions and their results include a K 
factor of 2. The purpose of this simple exercise is primarily to confirm that 
our numbers are not wildly divergent from mini-jet calculations. Our results 
are typically a factor of two larger (more at LHC) but this is easily understood 
because our results include all momentum modes. 

The number density of these out-of-equilibrium gluons can be related to 
the equilibrium entropy: S g i ue = 3.6-N g i ue . This is particularly so at the LHC, 
where because Q 2 S ^> Aq CD , elastic scattering dominates. The equilibrium 
entropy of gluons is, to within a factor of two, (which can be quantified in 
one's thermal+hydro model of choice), the entropy of pions. 
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